Low temperature transport through a quantum dot: 
finite-^/ results and scaling behavior 
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We calculate the conductance through a quantum dot weakly coupled to metallic leads, modeled 
by the spin-1/2 Anderson model with finite Coulomb repulsion U. We adopt the non-crossing 
approximation method in its finite-?/ extension (UNCA). Our results can be compared to those 
obtained with the exact numerical renormalization group method, and good agreement is found 
both in the high temperature (Coulomb blockade) and in the low temperature (Kondo) regime. 
We analyze the scaling properties of the low temperature conductance, and calculate the universal 
function which describes the electronic transport in the Kondo regime. Very good agreement with 
recent experimental results is found. Finally, we suggest a simple interpolating function which fits 
fairly well the calculated conductance in a broad temperature range. 



PACS: 72.15.Qm, 73.23.-b, 73.63.Kv 



I. INTRODUCTION 

The Kondo effect, a phenomenon discovered in the 
late 30s in diluted magnetic alloys, plays a crucial role 
in the low temperature properties of many strongly cor- 
related systems, such as Ce, Y and U compounds!!! Very 
recently Kondo-like phenomena were also observed in 
the low-tempjerature transport properties of quantum 
dot devices^u opening new opportunities to control the 
Kondo effect experimentally, and starting a new field of 
researchJiS 

A quantum dot device consists of a small sized quan- 
tum dot (QD) weakly connected by tunnel barriers to 
two electrodes, called source and drain. In this sys- 
tem the QD may be considered as an artificial atom, in 
which a well defined number of electrons, n, is confined. 
The energy required to add a new electron to the QD is 
U = e 2 /2C, where C is the capacitance of the QD itself. 
The energy U is determined by the Coulomb repulsion 
between two electrons in the QD and thus it scales with 
the inverse of the dot dimensions. Therefore, in small 
sized QD, U is usually larger than the coupling to the 
leads. An important aspect of a QD device is that the 
energy of electrons in the dot can be tuned by a gate 
voltage. Increasing the gate voltage, the energy neces- 
sary to add a new electron decreases, eventually making 
the addition possible. Thus it was shown that, at low 
temperature, the transport of electrons through the QD 
is allowed only at those values of the gate voltage at 
which the state with n electrons is realized and that with 
n + 1 electrons becomes-suddenly energetically accessible 
(Coulomb blockade). BE3 

The idea that Kondo-like phenomena should appear 
in such a system at very low temperature can be traced 
back to 1988. In that year it was recognized that the An- 



derson model, introc 



impurity in a meta 
pled to its leads .o 



jced in 1961 to describe a magnetic 
3 could also be applied to a QD cou- 
HE9 Soon afterwards several theoretical 
studies have been devoted to analyze the properties of 
such a system (see e.g. Refs. p~6] p3|) . It was predicted 
that at very low temperature (T <C Tk, where Tk is the 
Kondo temperature) a narrow peak should appear in the 
local density of states, close to the Fermi level. Thus, 
states belonging to opposite electrodes should mix easier 
than at high temperature (T 3> Tk), and the conduc- 
tance should increase. In addition, the Kaado-anomalics 
were predicted to appear only for odd and thus 

only if the total spin of the electrons in the QD, 5, is half 
integer. These predictions are now confirmed by several 
experimental resultsriQ The Kondo temperature of QD 
devices is however very small (usually less than a few 
hundred mK) compared to the one, of diluted magnetic 
alloys (usually a few K). Recently Jffl in some QD devices 
Kondo anomalies have been observed for even n as well; 
it is believed that these deviations from the even/odd n 
effect are related to the formation of integer spin states 
with 5 > 1. 

In the present paper we will focus our attention on 
those QD devices for which the Kondo effect occurs only 
for odd n. The transport properties of such a QD de- 
vice have been calculated by using many differeflt. ap- 
proaches, such as the equations of motion method,tZI the 
non-crossing approximation (NCA)lij'E3 and the 
numerical renormalization group (NRG) technique.LJLij 
Within the NRG it was possible to calculate the con- 
ductance as a function of the gate voltage, and thus to 
obtain - in the case 5=1/2- the even/odd effect in 
very good agreement with experimental results. Similar 
results were also obtained with the approximate equa- 
tions of motion method and with the NCA. The latter 
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was applied to low temperature transport of 
a QD in the infinite U limijt-and thus the even/odd alter- 
nation effect was obtainedE3 for n switching from n = 
to n = 1. 

Even if approximated, the NCA can be easily extended 
to models in which the realistic levels structure of the QD 
is kept into account (i.e. in the presence of several levels 
and possibly of an external magnetic field) , thus becom- 
ing a very convenient method for a direct comparison to 
transport properties of real QD devices. As a matter of 
fact, the NCA was used with success to study, e.g— cp-, 
bital degeneracy effects in diluted magnetic alloys.E-3c£l 
However, in the case of QD devices, it is first necessary 
to go beyond the U — ► oo condition, which does not 
account for the even/odd effect for all values of n. 

In the. present work we apply the finite- U NCA 
methodEZI to the QD device system, modeled by the spin- 
1/2 Anderson model. We calculate the conductance as a 
function of the position of the QD level, and study the 
system in the empty orbital, mixed valence and Kondo 
regimes. Our results are in good agreement with the ex- 
act numerical renormalization group ones. In particular, 
we obtain the correct behavior of the conductance both 
in the high temperature regime, T 3> 7jr, in which the 
conductance shows the Coulomb blockade peaks, and the 
low temperature regime, T < Tk, where the conductance 
increases for odd n due to the Kondo effect. 

The paper is organized as follows. In section II we 
describe the model Hamiltonian and give a short intro- 
duction to the finite-?/ NCA. In section III we show and 
discuss the numerical results, in particular the linear- 
response conductance, which we compare to exact results 
and to experimental data. In section IV we study the 
scaling behavior of the conductance as a function of tem- 
perature, and we suggest an empirical formula which is 
valid both in the Fermi-liquid and in the T ~ Tk temper- 
ature regime, where the conductance has a logarithmic 
behavior. 



II. MODEL AND METHOD 
A. Hamiltonian 



the QD and U is the Coulomb interaction among elec- 
trons in the dot; rida — d\d a is the number of electrons 
operator for a given spin in the QD; Vk a is the hybridiza- 
tion between the leads and the QD states, whose modulus 
is supposed to be k-independent, with |Vk ff | = Vgrm for 
(k, a) G S(D). The lead-dot coupling strengths are given 
by 

r s{D) =nV§ (D) J2^-^)- (2) 

k 

It was shownEl - through a unitary transformation of 
the band states - that the Hamiltonian (|l|) is equivalent 
to a two band Anderson model in which the first band 
does not interact with the QD, and the second is cou- 
pled to the quantum dot states through the hybridiza- 
tion V = y/Vg + Vp. Thus the problem of calculating 
the transport properties of Hamiltonian (pi) is reduced 
to the problem of calculating the spectral properties of 
the one band Anderson model, provided that the actual 
lead-dot coupling strength is given byc3 

r = r s +r D = 7rA(e F )^ 2 , (3) 

where N(ef) is the density of states (per spin) of the 
leads at the Fermi level, ep. It is reasonable to ap- 
proximate the conduction bands of the leads with those 
of a non interacting two dimensional Fermi gas. Thus 
the density of states (DOS) per spin may be written as 
N(e) — 1/2D, where D is one half of the band-width, 
therefore T = nV 2 /2D. 

The main difference between the Anderson model used 
in the ordinary Kondo problem and the one used here is 
the following. In the present model the energy difference 
£o — £f is not fixed, but, on the contrary, it can be tuned 
by a gate voltage V g , coupled to the QD through a ca- 
pacitor. From now on we set the Fermi level ef = 0. In 
first approximation, we can assume that — eo increases 
linearly with eV g (with common conventions for the sign 
of gate voltage). Thus the equilibrium thermodynamic 
properties of the system described by the Hamiltonian 
(|l|) are functions of two external parameters: the tem- 
perature T and the gate voltage V g . 



In order to describe a quantum dot (QD) coupled with 
its leads we adopt the Anderson Hamiltoniantj 

H = ^ ^kC^ffCktr + £Q nd-\n dl 

+ E (V^ctada + h.C.). (1) 

(k,<r)eS,D 

Here c^ CT (ck a ) creates (destroys) a conduction electron 
with momentum k and spin a in one of the two leads, 
which we label with S (source) and D (drain); S a (d a ) 
creates (destroys) an electron with spin a on the quan- 
tum dot; £q is the energy of a single electron localized on 



B. Linear response conductance 

In the linear response regime (Vsd V g , where Vsd 
is the source-drain voltage), the conductance of the sys- 
tem rQD-jf leads, G, may be written in a Landauer-like 
formB-ta 

G( ™ = x/_? (-M^ + i,)}) (-|) de, 

(4) 

where for simplicity we assume that the couplings to 
the leads are symmetric (Ts = I/d). Here / is the 
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Fermi distribution function and G R (e + irj) is the re- 
tarded local Green function, i.e. the Fourier transform 

ofG*(i) = -ie(t)({d(t),dt( )}> . 

The quantity Im {G R (e + i?7)} in Eq. (Q) is propor- 
tional to the local density of states, p(s), defined as 

p( e ) = -Ilm{G*( e + i77)}. (5) 

We point out that the effects of the Coulomb interaction 
among electrons in the QD are contained in the local 
density of states. The main purpose of this work is thus 
to calculate p(e) as a function of temperature and gate 
voltage, and the resulting linear conductance. 

C. Number of electrons 

The low temperature transport properties of a meso- 
scopic system such as the QD+leads previously described 
are characterized by the phenomenon of charge quanti- 
zation. When Coulomb blockade occurs the number of 
electrons in the QD is a fixed integer. This number can 
be changed by raising the voltage of the gate electrode. 
In this way the energy of the electrons in the QD is low- 
ered with respect to the Fermi energy in the leads. The 
change in energy necessary to add a new electron is ~ C/. 
This may be seen as the energy required to charge the 
QD, e 2 /2G, where C is the capacitance of the QD. Thus 
the average number of electrons on the QD, (n), with 

/+oo 
def(e)p(e), (6) 
-oo 

is an important quantity in this problem (here N is the 
degeneracy of £o an d in the present case N = 2). The 
number of electrons is — like the conductance — a func- 
tion of the local density of states. 

D. Finite-C/ non-crossing approximation 

In order to calculate p(e), the local DOS for the single 
impurity Anderson model, we adopt the fiaite-E/ non- 
crossing approximation (UNCA) approach,EZI an exten- 
sion of the non-crossing approximation (NCA) to the fi- 
nite U Anderson model. We point out that a finite U 
treatment is necessary in order to calculate the correct 
behavior of the conductance as a function of the gate 
voltage (and thus of the number of electrons on the QD). 

The NCA is a diagrammatic technique, intr oduced for 
the Anderson model in the U — >oo limit The main 
idea of the NCA is that the self energy can be expanded 
in a series of diagrams of order V 2 /N, where N is the 
degeneracy of the local level, usually large for Kondo im- 
purities in diluted magnetic alloys (e.g. N = 6 for Ce 
impurities). The non-crossing diagrams are summed up 
to all orders, and the first neglected diagrams are of cp-i 
der (V 2 /N) 2 . The NCA was proved to be successfully 



already for N — 2. Thus it was applied with good re- 
sults also to the,|Sein-l/2 U — ► oo Anderson model out 
of equilibrium JlaEJ e.g. it was used to calculate the be- 
havior of the conductance as a function of the chemical 
potential in the QD+leads problem. 

The NCA and the UNCA lead to a set of integral 
equations for the self energy, which have to be solved 
self-consistently. The number of equations depends on 
the number of many body states for the electrons in 
the QD. If U is finite there are four possible manybody 
states, namely |0, 0) (with energy e = 0), |1, J.) and |1, |) 
(degenerate, with energy e = eq), and |2,t|) (with en- 
ergy e = 2e + U). For each state a self energy and a 
ionic propagator are introduced. Thus we have three 
self energies, Yjq, Si, £2 and three propagators, Go, 
Gi, G2 (due to spin degeneracy Gij = G%i — Gi and 
£ij = Sij Xi). The UNCA equations may then be 
written asEZrEj 

2r f D 

Z (u> + ir)) = — de f(e)G 1 (cj + e + i V ) 7 (7) 

TT J-D 

r f D 

Ei(w + ir,) = - de [1 - /(e)]G (w -£ + «?) + 

TT J-D 

r r D 

+- / de f{e)G 2 {u J + e + ir 1 ) 1 (8) 

J-D 

or r D 

E 2 (cj + ir,) = — d£[l-/(e)]Gi(w-e + i»7), (9) 

TT J-D 

with 

G (w + ir)) = [lu + ir) - Z (lu + ir/)}- 1 , (10) 
+ iry) = [w + if] - s - Ei(w + ir?)] -1 , (11) 
G 2 {lo + iii) = [lo + irt - 2s - U - E 2 (^> + "7)]^ • (12) 

In the infinite U limit, the double occupation state is 
forbidden, and G2 may be neglected. 

Once this system of six equations is solved, it is possi- 
ble to express all the physical quantities in terms of the 
ionic resolvents. In particular the retarded local Green 
function may be evaluated by analytic continuation from 
the corresponding imaginary time propagator, which may 
be written asE3 

[G {z)G 1 (z + iu) + G 1 (z)G 2 (z + iu)] , (13) 
where Z is the QD partition function, i.e. 

Z = I e- z / k * T [Go (z) + G 1 (z) + G 2 (z)} (14) 
Jc 27ri 

Finally, the local density of states of the QD can be ob- 
tained from the retarded local Green function through 
Eq. @. 
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FIG. 1. Equilibrium density of states for the particle-hole 
symmetric Anderson model (U — — 2eo) at different temper- 
atures. Parameters: so = —1 meV, U — 2 meV, V = 0.35 
meV, D = 4 meV. For k B T = 2.6 x 1(T 4 meV (T = 3 mK) 



the local DOS at Fermi energy is nFp(0) — 0.'. 
the Fermi level ef is set equal to zero. 



In this paper 



The NCA and the UNCA break down below a tem- 
perature 2-,* ,-which in the Kondo regime is much smaller 
than .E-3~E3 Below this temperature spurious features 
show up, e.g., in the local density of states. It is known 
that the exact results-, ace recovered with the inclusion 
of vertex corrections Ji^nEl within a set of integral equa- 
tions which are numerically heavier than Eqs. ([7-12). 
This is, however, beyond the purpose of the present pa- 
per. Here we will use the UNCA without vertex dia- 
grams. We will see that, nevertheless, we can reproduce 
fairly well both the high and low temperature regimes, 
and that (above T*) the local DOS calculated with the 
UNCA is in good agreement with the exact numerical 
renormalization group (NRG) results. 



III. NUMERICAL RESULTS 

We solve the self-consistent UNCA equations, and cal- 
culate the local DOS, by using the fast Fourier transforms 
techniqueEj From p(e) we then obtain the Landauer con- 
ductance and the number of electrons as a function of the 
temperature and of the position of the QD level, Eq, pro- 



portional to —V g 



A. Local density of states 

The local density of states is shown in Fig. |l] as a 
function of the energy. The ratio — £q/T ~ 2.85 > 1, and 
thus the model describes the system in the Kondo regime. 
For simplicity, we show the results for the particle-hole 
symmetric Anderson model (U = —2eq). The Konde 
temperature may be estimated from Haldane formula," 



k B T K ~ (C/r/2) 1 / 2 exp [tc ( £o + U)/{2UT)] , (15) 

and in the present case we find k B Tx ~ 0.06 meV. 

In the high temperature limit (T ^> Tk, Fig. |l|a), the 
local density of states shows two broad resonances. The 
first one is located at energy e — E(n = 1)—E(n = 0)=£q, 
the energy required to put the first electron in the QD. 
The second resonance appears at e = E(n = 2) — E(n = 
1) = £o + U, the energy required to add an additional 
electron. The shape of the resonance peaks is Lorentzian. 
The origin of the peak broadening is mainly the coupling 
QD-leads measured by the width T, although a small 
thermal contribution is present since k B T>T. 

As the temperature decreases, a peak - which is the 
fingerprint of the Kondo effect - appears close to the 
Fermi level. The height of the peak increases on decreas- 
ing the temperature. Figure [j] shows that for T <C Tk 
the height of the Kondo peak tends to a maximum value, 
i.e. itTp(ef = 0) — ► 1, which is reached exactly only at 
T = 0. At k B T = 2.6 x 10~ 4 meV (T/T K ~ 0.01, Fig. 
0d) we find 7rrp(0) - 0.98. 

The evolution of p(e) as a function of temper ature , is, 
in fairly good agreement with the NRG results. Oc3£3 
In the very large U limit we find results consistent with 
those-obtained by using the NCA method with infinite 

u mm 



B. Linear-response conductance 



The linear response conductance G and the average 
numbeiujif electrons (n) are shown in Fig. |^ as a function 
of — £o-E3 For — Eq < U/2, when the number of electrons 
in the dot is n < 1, there are three relevant regimes of 
interest: the Kondo (K) regime T < —Eq < U/2, the 
mixed valence (MV) regime |eo| < T, and the empty 
orbital (EO) regime — £o < — T. In the opposite case 
— £o > U/2, when the number of electrons in the QD is 
n > 1, it is convenient to discuss the various regimes 
in terms of holes; thus the Kondo regime exists for 
U/2 < ~£q < U — r, the mixed valence regime for 
| — £q — U\ < r, and the empty orbital regime with n ~ 2 
occurs for — Eq > U + T. These regimes are schematically 
indicated in Fig. |^b. 

In the Kondo regime, we can extimate the Kondo 
temperature from Haldane formula, Eq. (|l5|). We find 
k B T K - 10~ 5 - 10" 2 meV for T - 0.094 (Fig. 
k B T K - 0.062 - 0.161 meV for T - 0.35 (Fig. '_ 

Figure ||a shows the results obtai^ied-for k B T ~ 0.5r ^> 
k B Tx- In this limit it was shownllinlZl that the conduc- 
tance has narrow peaks every time the average number of 
electrons in the QD increases by one. The peaks are sepa- 
rated by valleys in which the conductance is almost zexo. 
This behavior (which is due to the Coulomb blockadeO) 
is well reproduced in Fig. 0a. The peaks reach about 
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FIG. 2. Linear-response conductance G and average num- 
ber of electrons on the dot (n) , for two different set of parame- 
ters: (a) Parameters are: F = 0.094 meV, D = 4 meV, (7 = 3 
meV, and k B T = 4.3 x 10~ 2 meV (T = 500 mK) (Coulomb 
blockade); (b) Parameters are: F = 0.35 meV, D = 4 meV, 
U = 2 meV, and fc s T = 4.3x 10~ 4 meV (T = 5 mK). The var- 
ious regimes as a function of — eo are schematically indicated 
in (b) (see text) . 



e 2 /h, as expected from Coulomb blockade theory in the 
limit fc^T ~ r, and have a width of about 2I\ This 
broadening is due to tunneling, while the thermal broad- 
ening is negligible. The line shape is almost Lorentzian. 
If the temperature is raised so much that fc^T 3> F the 
height of the peaks becomes much smaller than e 2 /h and 
their width is controlled by thermal broadening only. 

Figure Hb^hows G and (n) for T -C Tk- In this regime, 
experiments3~l3 show that the valleys tend to raise when 
the number of electrons in the QD is odd, and remain 
almost unchanged when the number of electrons in the 
QD is even. In the very low temperature limit, the 
valleys .with an odd number of electrons evolve into a 
plateauQ at G ~ 2e 2 /h. This behavior is the manifes- 
tation of the Kondo effect in the mesoscopic transport 
properties of the QD+leads system. Figure ||b shows 
that when (n) ~ 1 the conductance calculated with the 
UNCA tends to G = 2e 2 /h, while a valley appears for 
(n) ~ and (n) ~ 2. Thus the Kondo regime is well 
described by the adopted method. The average number 
of electrons in Fig. ||b has a regular increase and does 
not have well defined plateaux as in Fig. [2^,, due to the 
larger value of T. 

Thus we have shown that the UNCA can describe well 




<n> = 2 



-2-10 1 



3 4 




-£ 







(meV) 



FIG. 3. Linear response conductance as a function of — eo 
(eo is always referred to the Fermi level, being ef — 0) the en- 
ergy of the dot level with respect to the Fermi energy, and for 
many temperatures, (a) Parameters: F = 0.35 meV, D = 4 
meV, U = 2 meV. Temperatures (from the bottom to the 
top curve): k B T = 0.43, 0.086, 0.043, 8.6 x I0" 3 , 4.3 x 10" 3 , 
8.6 x 10~ 4 , and 4.3 x 10~ 4 meV. (b) Parameters: F = 0.275 
meV, D — 4 meV, U — 2 meV. Temperatures (from the bot- 
tom to the top curve): k B T = 0.26, 0.17, 0.07, 0.035, 0.017 
8.6 x 10" 3 , 5.3 x 10" 3 , 2.6 x 10 
4.3 x 10~ 4 meV. 



3 , 1.3 x 10" 3 , 8.6 x 10" 4 , and 



both the Coulomb blockade and the Kondo effect. 
The results r-caa, be compared with the exact NRG 
calculations .EilH Agreement is very good in the Kondo 
regime, down to T ~ 0.01TV . In the mixed valence and 
empty orbital regimes the agreement decreases. At low 
temperature the UNCA breaks down and a spurious peak 
appears in the local DOS, close to the Fermi level, as dis- 
cussed for the usual NCA in Ref. Within the parame- 
ters and the temperature chosen in Fig. |^b, this happens 
approximately for — eo < 0.5 meV and — eo > 1-5 meV. In 
this region the conductance calculated with the UNCA is 
much larger than the exact NRG result. Therefore from 
now on we will plot the conductance only for the values 
of T and eo for which the UNCA is reliable. 

The linear-response conductance is shown in Fig. || 
as a function of — eo and for several different tempera- 
tures. The results are shown for two different choices 
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of the tunneling width: T = 0.35 meV (Fig. ||a) and 
r = 0.275 meV (Fig. ^b). In the Kondo regime, we esti- 
mate the Kondo temperature by using Haldane formula, 
and find k B T K ~ 0.062 - 0.161 meV for T = 0.35 meV 
and k B T K ~ 0.030 - 0.134 meV for T = 0.275 meV. 

Figure ||a shows that at high temperature (k B T > T ^> 
^bTk) there are two Coulomb blockade peaks at energy 
— So ~ and —£q ~ U, the maximum value of the conduc- 
tance being G ~ Q.5e 2 /h. As the temperature decreases, 
the Coulomb blockade peaks become closer to each other, 
and ideally merge into a plateau at T = 0. This is due 
to the fact that the energy — £q itself is renormalized by 
the many body effects, through the real part of the self 
energy. At very low temperature we see that G ~ 2e 2 /h 
(with maximum deviation of 10%) in the region in which 
the number of electrons is odd (n ~ 1). 

The same behavior is shown in Fig. ||b. Here we 
choose a smaller value of T, so that the Kondo temper- 
ature is smaller. At high temperature we see again the 
two Coulomb blockade peaks, located at — so ~ and 
— £o ~ U. At very low temperature the peaks approach 
each other and tend to the maximum value 2e 2 /h. How- 
ever in this case the deviation from the limiting value 
is larger (about 25%). This is due to the fact that, be- 
cause of the smaller Tk, down at T ~ 0.017V (about 
the lowest temperature that can be reached before the 
UNCA breaks down) we have 7rFp(0) ~ 0.75 instead of 
7rrp(0) = 1. Nevertheless, Fig. ||b reproduces the theo- 
retical NRG calculations in the Kondo regime at several 
different temperatures, see e.g. Fig . 3 in Ref. pi], Fig. 2 
in Ref. 



22, and Fig. 2a of Ref. E3l In addition, at high 
temperature Fig. |3Jb provides very good results in the 
whole range of — £q. 

It is worth to remind that when k B T > T the peak 
broadening is mainly due to thermal effects. Thus, when 
we start from the high temperature regime and lower 
the temperature, at first the peak width decreases. This 
effect is shown in Fig. On the contrary, when 

k B T ~ r the broadening.Js about 2T, in agreement 
with experimental results. au If the Kondo temperature 
is negligibly small, lowering the temperature below T/k B 
does not produce any further change in the peak width. 
However, if the Kondo temperature is finite, as soon 
as k B T k B Tx <C T the two peaks tend to become 
broader, and merge in a single large peak at zero tem- 
perature. This effect is shown in Fig. @. 

The UNCA results in Fig. | can be directly com- 
pared to the experimental measurements of the linear 
conductance vs the gate voltage. The curves in Fig. || 
reproduce very well the experiments described in Refs. 

in which the measurement of the Kondo effect in 
the equilibrium conductance of a QD+leads system is re- 
ported. Agreement is satisfactory for the temperature 
dependence and the broadening of the peaks, and it cov- 
ers the whole range of electron occupancy < n < 2. 
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FIG. 4. Universal behavior of linear-response conductance 
normalized to its saturation value Go- The conductance is 
displayed as a function of T/Tk, for three different values of 
£o/T chosen in the Kondo regime. Parameters: T = 0.275 
meV, D = 4 meV, U = 2 meV. 



IV. SCALING BEHAVIOR 
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In diluted magnetic alloys, the low temperature ( 
Tk) resistivity, p(T), follows a universal scaling law 
i.e. p/po = F(T/T K ), where p Q = p(T = 0) and F(x) 
is a function independent on system-related parameters. 
This happens because at low temperature and in the 
Kondo regime there is only one relevant energy scale, 
k B TK- Similar scaling properties have been repor 
the low temperature conductance of QD devices,! 
and carbon nanotubes.Ell For a QD device described 
by the spin-1/2 Anderson model the universal fpaction 
G/Go was recently calculated by using the NRGE3 Here 
we calculate the same function by using the UNCA, and 
compare our results with the NRG calculations and ex- 
perimental data. 



A. Universal curves 

The conductance calculated with the UNCA is shown 
in in Fig. | as a function of T/Tk- The different set 
of points correspond to different choices of — Eq/T in the 
Kondo regime (— £q/T = 3.6 corresponds to the particle- 
hole symmetric point, that is — £q = {7/2=1 meV in Fig. 
||b). The conductance is normalized to Go =G(T ,V g ), 
where To is the lowest temperature which can be reached 
before the UNCA breaks down. In order to achieve a uni- 
versal behavior we define the Konda ^temperature such 
that G(T K ) = Go/2, as often doneBE The relative de- 
viation from ideal result in the Kondo region for — £q/T, 
that is 5 = (Go — 2e 2 /h)/(2e 2 /h), depends on the choice 
of the parameters (we find 2% < \S\ < 25%). Figure | 
shows that the conductance follows the universal behav- 
ior in a fairly large range of temperatures. Departures 
from this universality are observed only for T 3> Tk , as 
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FIG. 5. Experimental data (Fig. 3c of Ref. 7) compared 
with the UNCA universal curve G/Go- Parameters of the 
theoretical curve: eo/T — -1-3.27, D = 4 meV, U — 2 meV. 
Here V g i is the gate voltage.0 




100 



FIG. 6. Universal curve G/Go (so/T = —3.6) compared 
with the interpolation function given in Eq. (^|) . The result- 
ing fit parameters are: a = l/ln(l + b + c) ~ 0.25, b = 43.4 
and c = 10.25. 



expected. We point out that at very low temperature 
(T/Tk <C 1) the-universal curve has the expected Fermi 
liquid behavior,L3 i.e. (G— Go) /Go oc T 2 , while at higher 
temperature (T/Tk ~ 1) the conductance is proportional 
to In (T/T K ). 

The universal curve of the conductance calculated with 
the UNCA is similar to NRG results for the ordinary 
Kondo problem (see e.g. Fig. 4 of Ref. || or Fig. 12 of 
Ref. |35|) and to NRG results for the conductance in QD 
devices (see Figs. 4 and 7a of Ref. p2h. 



B. Comparison with experiments 

The universal curve calculated here can be directly 
compared with existing measurements of the conduc- 
tance. In this section we will compare in particular with 
the experimental data taken from Ref. ^. In Fig. [| we 
show G/Go vs T/Tk, and plot both the experimental 
points and our UNCA curve. In Ref. ^ the experimental 
data were normalized to Go, the value of the conductance 
at the lowest temperature for which a measurement ex- 
ists, and the Kondo temperature was defined as the tem- 
perature such that G{Tk) — Go/2. 

We found the best agreements between UNCA and ex- 
periments for the choice £o/r = —3.27. Figure || shows 
that the crossover from the logarithmic behavior to the 
low temperature Fermi liquid regime is reproduced fairly 
well by the conductance calculated by using the UNCA. 



G(T) = G 



rpl 2 

1 K 



rp2 



rpl 2 

± K 



(16) 



where T' K = Tk 



2V* - 1 and G(T K ) = G /2. This for- 
mula reproduces well NRG results for T < Tk and has 
a single fitting parameter, s. However, in the interme- 
diate temperatures regime (T ~ Tk) the exact conduc- 
tance^ proportional to \h{T/Tk), as NRG calculations 
showJH3 and this behavior is not explicitly accounted for 
in formula (|l6|). Here we propose an alternative phe- 
nomenological formula, which (a) is still quite simple, 
(b) reproduces fairly well the calculated conductance in 
the whole temperature regime and (c) shows explicitly 
the logarithmic behavior in the intermediate temperature 
regime. Our formula has two fitting parameters instead 
of one. 

We notice that at very low temperature, in the Fermi 
liquid regime (T <C Tk), the conductance should show 
the Fermi liquid T 2 behavior, 



G = Go 1 



T 2 

1 K 



(17) 



where a is a parameter. Instead, at higher temperatures 
(T ~ Tk) a logarithmic behavior is expected 



Gcx Go In (T/Tk) 



(18) 



A function that satisfies both requirements is the follow- 
ing 



C. Interpolation formula 

In order to extract the Kondo temperaturje|Jpam ex- 
perimental data, it has become very commoroElO to fit 
the data with the empirical formula 



lr 1+oln 



(19) 



where a, b, and c are dimensionless parameters to be 
determined with best fit techniques. By defining Tk 
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as the temperature such that G(Tk) = Go/2, we find 
a = l/ln(l + b + c). Thus there are only two fitting 
parameters in our formula. 

In Fig. ^| we show the results of a fit for the choice 
— Sq/T = —3.6. The empirical formula reproduces ex- 
tremely well the UNCA results for the values b — 43.4, 
c = 10.25. These parameters slightly depend, of course, 
on the choice of — sq/T, and of the other UNCA param- 
eters. Thus Eq. ( pl^ ) reproduces the correct behavior, 
and it is physically correct both in the T -C Tk and 
T ~ Tk regime. We notice that the fourth order term is 
required to reproduce the behavior of the conductance in 
the intermediate temperature regime, while the T 2 term 
is required to reproduce the low temperature Fermi liquid 
behavior. 



V. CONCLUSION 

In the present work we have calculated the conduc- 
tance of a system made of a quantum dot coupled to two 
leads, described by the spin-1/2 Anderson model. We 
adopted the finite-t/ non-crossing approximation method 
(UNCA), which allowed us to calculate the conductance 
for the Anderson model with finite Coulomb repulsion. 
Thus we were able to study the conductance as a function 
of temperature and gate voltage. We have shown that 
the results obtained with this method are in good agree- 
ment with those obtained by ppips the exact numerical 
renormalization group (NRG).EirE3 In addition we repro- 
duced both the Coulomb blockade and the Kondo effect 
in quantum dots. Inclusion of a finite Coulomb correla- 
tion is important in order to describe correctly the exper- 
imental results in the whole regime of electron occupan- 
cies. The comparison with experimental data of Ref. [?] is 
fairly good, for what concerns both temperature and gate 
voltage dependence. Finally we have suggested a simple 
phenomenological formula which fits UNCA results both 
in the logarithmic and in the Fermi-liquid temperature 
regions, reproducing very well also the crossover between 
them. 

Although the spin-1/2 Anderson model can be solved 
by using the exact NRG, we believe that the UNCA is 
more suitable than NRG for extension to realistic sys- 
tems, e.g. for taking into account the effects of the elec- 
tronic structure of the dots, and thus we point out that 
the UNCA method can become an important tool to in- 
terpret the experiments. Possible applications are, for 
example, in the explanation of the Kondo effect in quan- 
tum dots fownteger spinEI and the Kondo effect in carbon 
nanotubesO 
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